Transition Radiation of Moving Abrikosov Vortices. 



O 

o 
o 

(N 

in 



o 
o 

Oh' 



^3 

o 

> 

in 

(N 
(N 

o 

oo 
On 
-i— > 



i 

c 

o 
o 



X 



O.V. Dolgov 
Max- Planck-Institut fur Festkorperforschung, 
Heisenbergstrasse 1, D-70569 Stuttgart, Germany 



N. Schopohl 

Eberhard-Karls-Universitat Tubingen, Institut fur Theoretische Physik, 
Auf der Morgenstelle 14, 
D-72076 Tubingen, Germany 

We predict that Abrikosov vortices moving towards the surface of a superconductor emit electro- 
magnetic radiation into free space. Passage of a single vortex line travelling at constant speed v 
from a superconducting half space, x < 0, towards the vacuum half space, x > 0, is discussed in 
detail. The frequency distribution of the radiated intensity displays a pronounced maximum at 
microwave frequencies around uj ~ \v x \/X, where A is the magnetic penetration depth. Coherent 
motion of a lattice of flux lines leads to constructive interference and increases the strength of the 
radiated power per vortex line by a large factor. Also the radiation source may be used to measure 
the velocity of Abrikosov vortices and to determine the magnetic penetration depth. 
PACS numbers: 41.20.Bt, 41.60.-m, 74.25.Nf, 74.60.-w 



A charged particle, moving inside a homogenous 
medium at constant speed v and not exceeding the phase 
velocity c of light in the medium, cannot radiate elec- 
tromagnetic waves. However, when the particle passes 
from one medium into another with different electromaa,- 
netic properties there may occur transition radiationtl,0. 
In sharp contrast to Cherenkov radiation there exists no 
lower bound on the velocity of a charged particle for tran- 
sition radiation. 

In the following we suggest that the phenomenon of 
transition radiation should also occur for electromagnetic 
objects like Abrikosov vortices moving towards the sur- 
face of a superconductor. In our calculation we assume 
all vortex lines being orientated parallel to the z-axis and 
we assume z • v = 0. rj-. 

Let us first consider a single Abrikosov vortex at resfffl. 
Such a vortex is a topological object inside an ordered 
medium, namely a line defect in a superconductor. Also 
it is a quantum object. The gradient of the phase of 
the pairing wave function drives Cooper pairs ( i.e. pairs 
of correlated electrons near the Fermi surface carrying a 
small center of mass momentum) in circles around the 
vortex axis, while the modulus of the pairing potential 
vanishes at the center of the vortex core (say the node is 
at position R in the xy-plane) . The flow of mass around 
the vortex axis displays a quantized circulation in units 
of tt" - , and the magnetic flux associated with such a line 
defect is quantized in units of $ = the quantum 
of flux. Due to Meissner screening the magnetic field 
around the vortex decays exponentially fast on the scale 
of the London penetration depth A with increasing dis- 
tance |r — R| 3> A to the vortex core: 
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In sharp contrast to this fast decay of the magnetic field 
B(r) is the asymptotic behaviour of the vector potential 



A(r) around an Abrikosov vortex, which is connected to 
the magnetic field via B(r) = rot A(r): 



A(r) = ^zA 



R 



|r-R|" 



for |r - R| > A 



(2) 



Note that A(r) is independent on A far outside the 
London core (apart from exponentially small correction 
terms). 

This slow decay of the vector potential around a sin- 
gle Abrikosov vortex is reminiscent to the behaviour 
of the vector potential in the exterior region of a long 
solenoid, like in an Aharonov-Bohm experiment. Con- 
sider a solenoid consisting out of a piece of infinitesi- 
mally thin wire wound helically around the axis z of an 
infinitely long cylinder with radius r a . Through the wire 
of the solenoid flows a suitable electrical current I a to 
generate a constant magnetic field B^ = B^z inside, 
but outside the solenoid the magnetic field is zero. The 
vector potential A( CT ) generated by the solenoid in its ex- 
terior region is the gradient of a 'phase': 



A^(r) = l 4^ A (r-R) 
U 1 |,W(r-R) 



for 
for 



|r - R| < r, 
Ir - Rl > r„ 



where e^ = f±M , W(r) = zA 



(3) 

and r = (x, y). 



Consider a closed path T in the exterior region of the 
solenoid with winding number 1 around the z— axis. In- 
tegration of A( CT )(r) along T leads to: 

<fdr- A (CT) (r) = $ CT = / d 2 r 2 • rot A (<t) (r) = nrl B {a) 

(4) 

(apply Stokes theorem to the area enclosed by T) . So, 
the constant Q a is equal to the flux passing through the 
solenoid. Adjusting the control current la- such that $ CT 
becomes equal to one quantum of flux, | <E> CT | = <&o , and 
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taking the limit r a — > + , the thin solenoid produces 
(in its rest frame) a singular electromagnetic field at its 
center: 



B^(r) = $ a z6 (2 \r R) 
EW(r)=0 



(5) 



Consider now a thin solenoid travelling relative to the 
lab frame at a constant speed v along the trajectory 
R = vt. Such a solenoid is surrounded by a time de- 
pendent Aharonov-Bohm type vector potential giving rise 
not only to a singular magnetic field B^^r, t) but also 
to a singular electric field EW(r,(): 



B^(r,t) = $ CT z^ 2 )(r- vt) 
E< ff >(r,t) = --ABW(r,() 



(6) 



It is consistent (and convenient) to associate with the 
moving thin solenoid a conserved charge-current density 
via: 

47rp (CT) (r,<) = divE (<T) (r,t) (7) 

— J (<T) (r,t) = rotB^(r,f) - -d t E,^(r,t) 

c c 

Since in the geometry under consideration translational 
invariance is broken at x — it is natural to require that 
no current must flow from the region x < (supercon- 
ductor) into the region x > (vacuum): 



limn- J (ff) (r,t) = 



(8) 



(n denotes the normal unit vector of the interface). By a 
simple trick we are able to fulfill this boundary condition 
(which actually avoids solving an integral equation prob- 
lem of the Wiener-Hopf type). Let us superpose to the 
field around a single flux line situated at its (instanta- 
neous) position R = vt another singular field generated 
by a mirror vortex (with opposite magnetic moment) sit- 
uated at the mirror reflected position R = vt , where 
v = (—v x ,v y ). The resulting fields (inside the region 
x < 0) are then: 



B (<j) (r, t) = $ CT z [<5 (2) (r - vt) - <5 (2) (r - vt) 



(9) 



E ((T) (r,t) = $ ff 2A 



I*< 2 >(r-vt)--*C2)( r _vi) 



Inserting these expressions for B' CT ) and F,^ into Eq.(Q) 
there results a conserved charge-current distribution 
which generates the asymptotics of the Aharonov-Bohm 
type vector potential around a moving Abrikosov vortex 
line far outside the London core region. By construction 
this charge-current distribution also fulfills the boundary 
value condition Eq.(||). 

Let us next discuss the effect of Meissner screening. 
Since a superconductor is an ordered medium it will re- 
act to a specified external influence by a specific response 



that takes into account the formation of induced currents 
and of induced polarization charges inside the sample. 
This means a parametrically driven flux line moving at 
speed v inside a medium with Meissner screening (imag- 
ine an external agent who pulls the solenoid along) is 
surrounded by total current and charge densities of the 
form 



j(* ot ) (r, u) = J (md) (r, u) + J (ff) (r, u) (10) 



^ (r, w) - p {md) (r, w) + p {a) (r, u) 



Here j( md ) and represent the induced response of 

the ordered medium to the applied drive. Explicit ex- 
pressions for the drive terms (r, to) and ^ ff '(r,u) 
follow inserting Eq.(||) into Eq.(^) and calculating the 
Fourier transforms with respect to time. Since the in- 
duced response is retarded in time it is convenient to 
switch to the frequency representation. Suitable expres- 
sions for p( m<i ) (r, u>) and j( md ) (r, uj) follow from standard 
linear response theory. Since any local perturbation in- 
side the bulk of a superconductor heals typically on the 
length scale of the coherence length £ of the Cooper pairs 
(£ <C A in extreme type-II materials), and since the phe- 
nomenon we study is related to low frequencies in the 
microwave regime, a suitable (local) expression for the 
induced charge-current distribution is: 

j(mrf)(r )CJ ) = _JL E^fr.w) (11) 

47r iu A, , 



p( md )(r,w) = — divJ ( " ld) (r,w) 



In a Hwo fluid' type description of the electromag- 
netic response of a superconductor it is found that the 
length A w is a frequency dependent (but local) penetra- 
tion depth, that interpolates between the skin penetra- 
tion depth, which is characteristic for the normal metallic 
state at frequencies above 2 A, to the magnetic penetra- 
tion depth (MPD) at zero frequency: 



^_^e(M-2A) + -Ie(2A-M) (12) 
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Here u> p i — \J ilT J^* 2 is the bulk plasma frequency, 7 de- 
scribes quasiparticlertdamping, A = A (7) denotes MPD 
of a superconductor^. In the clean limit (and at zero 
temperature) lim 7 ^ A (7) = ^- = A/,, the standard def- 
inition of the London penetration depth. Of course, the 
nonlocality of the electromagnetic response kernel be- 
comes eventually important very near to the surface (e.g. 
anomalous skin effect). 

In the presence of Meissner screening the (total) elec- 
tromagnetic field around a travelling thin solenoid con- 
sists out of two parts. The first part is entirely due to 
the Meissner screening currents which flow in the exte- 
rior region of the solenoid. This part actually represents 
(within the accuracy of local London electrodynamics) 



2 



the electromagnetic fields around a moving Abrikosov 
vortex, provided the flux inside the solenoid is fixed with 
the opposite sign: <fr a = — <3>o- The second part origi- 
nates from the infinitesimal small interior region of the 
thin solenoid, where there exists (in the rest frame of 
the solenoid) an artificially fixed singular magnetic field: 
£j(<r) _ _±<t Thjg feature of the thin solenoid model is 

in sharp contrast to what really happens at the center 
of an Abrikosov vortex, where the pairing wave function 
displays a node and where the magnetic field assumes a 
finite and (almost) constant value inside the inner core of 
the vortex , a region of size £ -C A. As a result, the self 
fields E^) (r, uj) and B^) (r, uj), which originate from the 
interior region of the solenoid, should be subtracted from 
the total fields E(* ot) (r,w) and B( tot )(r,a;) in order to 
determine the physical fields E(r, uj) and B(r, uj) around 
an Abrikosov vortex line: 



B( to ')(r,w)= B(r,w) +BW(r,u) 
E( to ')(r,w)= E(r,u)+EW(r,u) 



(13) 



The fields ^ tot) (r, uj) and B( tot > (r, uj) are solutions to the 
Maxwell equations in the presence of the total charge- 
current distribution specified in Eq.(|10|): 



47T 

c 



3 {a) (r,uj) + J( md )(r,w) 
= rotB( tot )(r, W ) + ^E( tot )(r, W ) 



(14) 



Upon insertion of Eq. (JTl]) this transforms into: 

4tt _/-„■,, , c 1 
iuj A 2 



c 
rot 



J (,T) (r,w) 
B(r,cj)+B (CT) (r,w) 



E(r,w) + ~E {a) {r,uj) (15) 



iuj 
c 



E(r,w) + E (ff) (r,w) 



Next we insert J^) from the definition Eq.(f7f) and obtain: 



second equation reduces to the standard London equation 
determining the. magnetic field around an Abrikosov vor- 
tex line at resta,0ja. Note that for simplicity we assume 
here that there exists no external static magnetic field in 
vacuum, i.e. lim^^o B(r,w) = for x > 0. 

Next we discuss scattering solutions to our dynami- 
cal London equations, which describe the transition ra- 
diation effect for moving Abrikosov vortices. Before we 
expand details of the new radiation source let us com- 
ment on an. interesting difference to standard transition 
radiationn,a. In the case of standard transition radiation 
the moving particle, say an electron or a proton, actu- 
ally exists on either side of the interface passing from 
one medium into another with different electromagnetic 
properties. In the case of a moving vortex the 'particle' 
(i.e. the travelling singular line representing the center of 
the vortex), exists only inside the superconductor, but 
ceases to exist on the other side of the interface (in vac- 
uum). 

A Fourier transformation of the dynamical London 
equations with respect to the y— coordinate simplifies the 
problem and leads to coupled ordinary differential equa- 
tions. It is sufficient to consider the tangential compo- 
nent E y {x) — Ey(x,k y ,Lu) of the electric field, since the 
normal component E x (x) = E x (x, k y ,ui) follows from 



knowledge of d x E y 



From both components of the 



electric field the magnetic field B z (x) — B z (x, k y ,u>) may 
be calculated using Faraday's law. Guided by the struc- 
ture of our problem we make now the following ansatz 
for E y (x) in terms of three unknown coefficients as , bs 
and ay , and in terms of three (complex) wave numbers 
qx , qoo and K u : 



E y {x) 



a s exp (q x x) + b s \e° 

ay exp (—qoox) 



-iK u x\ 



(18) 

for x < 
for x > 
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iuj A 2 



E(r,w) 



E 



(r, uj) = rot B(r, uj) H E(r, uj) 

c 

(16) 



Using the Faraday law 
rotE(r,w) 



— B (t,uj) 
c 







it is not difficult to derive dynamical London equations 
for the physical electromagnetic fields around a moving 
vortex line: 



rotrotE(r,w) + E(r,w) 



1 

A~ 
1 

A? 



^E('>M (17) 
B(r,w)= ^B^M 



For x > we have p- = 0, of course, since there is noth- 
ing in vacuum. In the static limit, |v| — * and uj — > 0, the 



<1\ 



1 1' 2 

K y 



l 

A^ 



K„ 



kyVy 



(19) 



At infinity (|x| — > oo) we pose the usual radiation bound- 
ary conditions. So, in order to avoid exponential growth 
for x — ► ±oo , the branch cut of the square root is de- 
fined such that Re q\ > and Re goo > 0. The coefficient 
bs follows from the requirement that E y {x) should be 
a solution to the differential equation in the half space 
x < 0. The coefficients as and ay follow from the stan- 
dard boundary value conditions for the tangential and 
normal components of the electric field at x = 0. The 
results are: 



sign(Uz 



A 2 , 



{^ qx ) 2 -^ ky K u 
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( l'x 
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as 



ay 



2b s %r qx 



2b s 



(20) 



The coefficient 6s displays a resonance around uj = lo p i, 
while the coefficients as and ay display also a resonance 
around uj = the surface plasma frequency. The so- 
lution for Ey(x) in vacuum (x > 0) describes an out- 
ward travelling wave front, representing electromagnetic 
radiation. The radiation part is entirely determined by 
oscillating terms, which exist only in the range of wave 

numbers \k 



< 



, so that q 



-£sign(u;)y^ 



/■2 



is purely imaginary. For \k y \ > |^| the modulus of the 
electromagnetic field in vacuum is exponentially small for 
x » -rjjT , and may be safely neglected. 

Consider now a planar detector screen at a distance 
in the far zone in vacuum, with the plane 



x T > 



M 



of the screen orientated parallel to the superconductor- 
vacuum interface (at x = 0). The rate at which the 
emitted radiation is incident onto a target area element 
n da — —x da at position yt — (xt,Vt) on the screen 
may be calculated from the associated Poynting flux 
j- [E {tt 1 t) A B {tt, £)]•(— x). Let us denote the amount 
of energy c?£(r-r) so radiated (per length unit L of the 
vortex line) onto an infinitesimal target area element 

dax = L ■ dyx by de }^ ■ Then the total energy e ra d 
radiated into free space (per unit length L of the vortex 
line) is obtained integrating ^ over time and over the 
the entire screen area (wide detector): 



£rad 



c 
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duj 
2^ 



dk 

— — E v (x T , k y ,Lo)B z (x T , 

Z7T 



Let us introduce scaled variables: uj = 



Ax, 



ky — \hk y , q 



= Al^oo , q\ = Al^a and v x = v y — 
— , so that energy is measured in units of fku p i , length is 
measured in units of A^ = and velocity is measured 
in units of c. We may then rewrite the vacuum amplitude 
ay = a (k y ,uj) = Ay (k y ,Zj) in dimensionless form: 



A 



V - =2 



sign(v x ) 



ql 



K y 



'JykyU 



K {U - Vykyf + V 2 x q\ 



d 1 I q c 



u 2 q\ 



It follows 



2ttA, 



i n de 

2./_ a T., (1 



(21) 



doj uj \ Ay (uj sin 8, uj)] 2 



(22) 



The function I (8 : U) — uj ■ \ Ay (uj sin 8 : uJ)\ is the 
angle-frequency distribution of transition radiation emit- 
ted into free space by a single moving Abrikosov vortex. 
Apparently, the dependence of I(8,uj) on the angle 9 
for a non relativistic speed is such that it is practically 
constant for all directions, with the exception of grazing 
angles around \6\ ig ^. In Fig.(l) we plot the radiated 

intensity I(uj) = f\ ^-1(8, ZD) vs. scaled frequency uj 

for various vortex speeds. Values for the gap, the quasi- 
particle damping 7 and the vortex speeds are taken as 
indicated in the Figure caption. The radiated intensity 
I(ul) displays a sharp peak at a frequency 



2nu 



1 



'peak 



^peak 



V3A( 7 ) 



(23) 



This is in agreement with the physical picture, that it 
takes a time T ~ for the London core of a vortex 
to pass the interfacial region. Our analytical estimata- 
tion for uj pea k applies provided the frequency of the emit- 
ted radiation fulfills Hcu < 2A. Indeed, even for a very 
high vortex speed around \v x \ ~ Vp on finds hu pea k = 
1 '" •'• - *sea • A < A , since £ RC s - % ■ A/.. 



Rather high vortex speeds around v x ~ lO 3 ^ have been 
measured in layered-type-II materials in the experiments 
of Doettinger et alE. For a London penetration depth 
of order A (7) ~ 10 3 A we expect the position of the 

peak at v peak = ~ ~ 1 GHz , a typ- 

ical microwave frequency. Also the height of the peak 
intensity increases linearly with the vortex speed v x , as 
shown in Fig.(l), in agreement with our analytical esti- 
mate I {uOpeak) — j\/3^-. Finally, the total amount of 
energy that gets emitted into free space due to the tran- 
sition radiation effect follows integrating over all frequen- 
cies. Certainly, at very high frequencies, say for uj > to p i , 
our local electrodynamics becomes invalid. But since the 
integral converges rapidly the error is small if we simply 
extend the upper integration limit towards infinity. We 
obtain: 



E-rad 



$0 



2ttA (7) 



(24) 



Assuming Al ~ 100 nm and a vortex speed of order 
'^f ~ 10~ 5 gives e rad ~ 10 -22 W ■ sec (per length 
L = 1 mi) . Due to the large impedance mismatch of 
the superconductor with the vacuum (Zq ~ 377 Q) it is 
not sumcising that only a tiny fraction of the nucleation 
energyEl,Q of a single Abrikosov vortex gets radiated into 
free space. So, when the line is destroyed during its pas- 
sage through the interfacial region, the energy that is not 
radiated into free space gets converted into elementary 
excitations (e.g. polaritons) and, eventually, into heat. 

A larger effect arises when several vortices, say a bun- 
dle of altogether |A| flux lines, move coherently at con- 
stant speed v towards the surface of a superconductor. 
According to the superposition principle the (effective) 
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drive term for several flux lines, situated at time t = 
at positions yl — {xl,Vl) inside the superconductor, 
follows from Eq.@ upon summation over the instanta- 
neous positions R^ = yl + v t of the individual vortices 
(and the instantaneous positions of the mirror vortices at 
Hl = tl + v t) . The linearity of the dynamical London 
equations Eq.(l^) implies that the electromagnetic field 
E(A) generated by such a bundle of flux lines is: 



E 

r L eA 



e -H.kyVL+K,,*L)E(x, ky,0j) (25) 



From the theory of diffraction of electromagnetic waves 
it follows that the associated Poynting flux, and hence 
the radiated peak intensity emitted by a periodic ar- 
rangement of coherently moving flux lines (for example, 
a square vortex lattice A, with lattice constant d\, con- 
sisting out of |A| = N x ■ N y flux lines), should be in- 
creased by a factor proportional to the number of partic- 
ipating modes compared to the case where only unlocked 
modes are superposed. The total radiated energy (per 
unit length L of the vortex lines) is 
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'rad 
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The structure factor of the moving vortex lattice 



S A (k y ,uj) 



r L eA 



(26) 



represents, for large N x and N y , a series of delta spikes, 
selecting wave numbers k y and frequencies u> such that 
k y = ^ny and cj = (n x v x + n y v y ), with n x and n y in- 
teger. For intermediate magnetic fields B c \ B c2 
the lattice constant g?a is much larger than the coherence 
length £, but smaller than the London depth Al , so that 
£ <C d\ <C Xl- We derive from this a simple estimation 
for the radiation energy emitted into free space by alto- 
gether |A| 3> 1 coherently moving vortex lines (arranged 
in the geometry of a square lattice as described above): 



_(A) 

'rad 
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7T 



d\ 



\Vx_ 

c 



(27) 



The principal result is that the transition radiation ef- 
fect for a cluster of altogether | A| 3> 1 coherently moving 
flux lines is a linear function of the vortex speed \v x \. The 
emitted power per vortex is huge compared to the radia- 
tion emitted by a single vortex by a factor proportional 
to the number of vortices, since -4- oc |A|. This enhance- 
ment is due to constructive interference, like in the super 
radiation effect of laser physics. Assuming Xl — 100 nm 
and g?a = 0.5 Xl gives a total of |A| ~ 0.25 x 10 12 vor- 
tices per cm 2 . This leads (per length L — lem) and 



a vortex speed v x 



10" 



to a value of or- 



der £ 



(A) 



that a sensitive S-I-S Josephson radiometer could be used 
to detect the predicted spectrum of transition radiation 
in free space. The sharp peak of the radiated intensity 
around u ~ allows to measure MPD for a known 
vortex speed, and vice versa. 

Our calculations indicate that transition radiation is 
a generic phenomenon, that should exist in other types 
of ordered condensed matter containing mobile (topologi- 
cal) defects. We expect transition radiation to be emitted 
into free space by rapidly moving dislocations, by mag- 
netic bubbles and Bloch walls in magnets, and also by 
travelling polarons. 
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FIG. 1. Frequency distribution 7(ZZ?) of intensity of tran- 
sition radiation emitted into free space for different vortex 
speeds; the labels a),b),c),d) correspond to v y — and to 
v a /c = 2 x 10~ 5 ,1.5 x 10~ 5 ,10~ 5 ,5 x 10~ 6 , respectively. 
Assumed values for scattering rate j/uipi ~ 10 -4 and gap 
A/u; pl = 10" 3 . 



0.12 x 10 i W sec. This estimation suggests 
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